Coupled cluster calculations of ground and excited states of nuclei 
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The standard and renormalized coupled cluster methods with singles, doubles, and noniterative 
triples and their generalizations to excited states, based on the equation of motion coupled cluster 
approach, are applied to the 4 He and ls O nuclei. A comparison of coupled cluster results with the 
results of the exact diagonalization of the Hamiltonian in the same model space shows that the quan- 
tum chemistry inspired coupled cluster approximations provide an excellent description of ground 
and excited states of nuclei. The bulk of the correlation effects is obtained at the coupled cluster 
singles and doubles level. Triples, treated noniteratively, provide the virtually exact description. 



The description of finite nuclei requires an understand- 
ing of both ground- and excited-state properties based 
on a given nuclear Hamiltonian. While much progress 
has been made in employing the Green's Function Monte 
Carlo 1] and no-core shell model @ techniques, these 
methods have apparent limitations to light nuclei. Given 
that present nuclear structure research facilities and the 
proposed Rare Isotope Accelerator will continue to open 
significant territory into regions of medium-mass and 
heavier nuclei, it becomes imperative to investigate meth- 
ods that will allow for a description of medium-mass sys- 
tems. Coupled cluster theory is a particularly promising 
candidate for such an endeavor due to its enormous suc- 
cess in quantum chemistry 0, S IE S S • 

Coupled cluster theory originated in nuclear physics 
[lfl lllj around 1960. Early studies in the seventies 0] 
probed ground-state properties in limited spaces with 
free nucleon-nucleon interactions available at the time. 
The subject was revisited only recently by Guardiola et 
al. 0|, for further theoretical development, and by Mi- 
haila and Heisenberg 0] , for coupled cluster calculations 
using realistic two- and three-nucleon bare interactions 
and expansions in the inverse particle-hole energy spac- 
ings. However, much of the impressive development in 
coupled cluster theory made in quantum chemistry in 
the last 15-20 years [EI El HE El still awaits applications 
to the nuclear many-body problem. 

In this Letter, we apply quantum chemistry inspired 

coupled cluster methods |lllHll0,|EHH[lB] t0 finit ° 
nuclei. We show that the coupled cluster approach is 
numerically inexpensive and accurate by comparing our 
results for 4 He with results from exact diagonalization in 
a model space consisting of four major oscillator shells. 
For the first time, we apply coupled cluster theory to 
excited states in nuclei, exploiting the equation of motion 
coupled cluster formalism 0, [h|. We discuss several 
approximations within coupled cluster theory and also 
compute the ground state of the 16 O nucleus within the 
same model space. We remind the reader that certain 



acronyms have become standard in quantum chemistry. 
For this reason, we use the same abbreviations in this 
Letter. 

Coupled cluster theory is based on an exponen- 
tial ansatz for the ground-state wave function |^o) = 
exp(T)|$). Here T is the cluster operator and |<&) 
is the reference determinant. In the CCSD ("coupled 
cluster with singles and doubles") method, we trun- 
cate the many-body expansion of the cluster opera- 
tor T at two-body components. The truncated clus- 
ter operator T (CCSD\ used in the CCSD calculations, 
has the form 17): T( CCSD ) = T x + T 2 . Here T x = 

Ei, 4 a ° a * and T 2 = lI2 lJ ,ab t ab aaaba j a * are the singly 
and doubly excited clusters, with indices k (a,b,c) 
designating the single-particle states occupied (unoc- 
cupied) in the reference Slater determinant |$) and 
a p (a p ) representing the creation (annihilation) oper- 
ators. We determine the singly and doubly excited 
cluster amplitudes t\ and t^ b , defining T\ and T 2 , re- 
spectively, by solving the nonlinear system of algebraic 
equations, ($<>|i?( CCSD )|$) = 0, ($^|i?( CCSD ) |$) = 0, 
where ij(ccsD) = ex p(-TC CCSD )) exp(T( CCSD )) is the 
similarity-transformed Hamiltonian and |$") and |$° & ) 
are the singly and doubly excited Slater determinants, 
respectively. Once T\ and T 2 amplitudes are determined, 
we calculate the ground-state CCSD energy Eq CCSD ^ as 
($|fl"( CCSD )|$). For the excited states \^k) and energies 
^(ccsd) ^ K > Q ^ we apply the EOMCCSD ("equation 
of motion CCSD" ) approximation, in which 



**> =i?£ CSD) exp(T( CCSD >)|<&). 



(1) 



Here I? 



(CCSD) _ 
K ~ 



Rn 



Ri + i?2 is a sum of the refer- 



ence (Rq), one-body (Ri), and two-body (R 2 ) compo- 
nents that are obtained by diagonalizing the similarity- 
transformed Hamiltonian ff( CCSD ) in the same space of 
singly and doubly excited determinants |<5") and HJ^l 6 ) 
as used in the ground-state CCSD calculations (lilllq. 
The CCSD and EOMCCSD methods are expected to 
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describe the bulk of the correlation effects with inex- 
pensive computational steps that scale as n„n„, where 
n (n u ) is the number of occupied (unoccupied) single- 
particle orbitals. While the inclusion of triply excited 
clusters T3 and three-body excitation operators R3 in- 
creases the accuracy of the method, the resulting full 
CCSDT ("T" stands for "triples") [Hi and EOMCCSDT 
[l9j | methods scale as n^n^ and are rather expensive. For 
this reason, we add the a posteriori corrections due to 
triples to the CCSD/EOMCCSD energies, which require 
n^n 4 noniterative steps. The ground- and excited-state 
triples corrections, Sq and 8k (K > 0), respectively, are 
calculated with the CR-CCSD(T) ("comp letely renor- 
malized CCSD(T)") approach 0, H E3, Ef in which 

S i< = m E ^K\^)Mf c (K)/A K (K>0). (2) 

ijk,abc 

Here are the triply excited determinants and 

■M^ bc (K) are the generalized moments of the CCSD 
(K = 0) and EOMCCSD (K > 0) equations [U HJ |H , 



M 



^w = « c i^ (CCSD) 4 CCSD) i^ 



(3) 



where S, 



(CCSD) 





= 1 and S 



(CCSD) _ ^(CCSD) 



K 



K 



for K > 0. 



They can be calculated using the CCSD and EOMCCSD 
cluster and excitation operators T^ CCSD ^ and R^ C L> 



respectively. The Ak denominators are defined as 



A A - = (*a-|4 CCSD) exp(T 



n(CCSD) 



1$) 



(4) 



The states \&k) m Eqs. © and 10} include the lead- 
ing triples contributions resulting from the perturba- 
tive analysis of the CCSDT and EOMCCSDT equa- 
tions. We have | * ) = Pexp(T( CCSD ) + T 3 )|$) and 
\9 K ) = P(R^ CSD) + J R 3 )exp(T( CCSD ))|$) for K > 0, 
where P is a projection operator on the subspace spanned 
by the reference $) and singly, doubly, and triply excited 
determinants. The most complete forms of T3 and R3 
denning the CR-CCSD(T),c approximation are HEU 

^ = m E (Mf c (0)/Dt%(0))a a a b a c a k a jai , (5) 

ijk,abc 

^ = M E (MT c (K)/D^(K))a a a b a c a k a jai ,(6) 

ijk,abc 



where Df^(K) 



^ CCSD) -(^I^ (CCSD) I^)-Inthe 
case of the ground-state calculations, we also consider 
simplified variants of the CR-CCSD(T) theory, termed 
CR-CCSD(T),a and CR-CCSD(T),b. In the case of CR- 
CCSD(T),b, the perturbative denominator D?M0) is re- 



placed by 



(foabcl rr(CCSD) 
\ ijk I 1 



$^ fe c ), where H\ 



'ijk < 
(CCSD) 



is the 



one-body part of #( CCSD ). For CR-CCSD(T),a we re- 
place £>"ife(0) by the standard many-body perturbation 



theory (MBPT) triples denominator (e-i + 



£b — e c ), where ti and e a are the diagonal elements of 
the Fock matrix. Very accurate results for the excitation 
energies Ek — E of many-electron systems are obtained 
if we use the complete CR-CCSD(T),c theory to calcu- 
late the energies of excited states and the CR-CCSD(T),b 
approximation for the ground-state energy [2l|. For the 
ground states, it may sometimes be worthwhile to re- 
place the Ao denominator, Eq. which rcnormalizcs 
the triples correction <5o by 1, since Ao equals 1 plus terms 
of the second MBPT order or higher [20j. We indicate 
this by using acronyms, such as CR-CCSD(T),c/A = 1 
(as opposed to CR-CCSD(T),c, where A is included) 

We use the Idaho-A nucleon-nucleon potential 
which was produced using techniques of chiral effec- 
tive field theory [2~1. l25j . Modern two-nucleon inter- 
actions, such as Idaho-A, include short-range repulsive 
cores that require calculations in extremely large model 
spaces to reach converged results In order to re- 

move the hard-core part of the interaction from the prob- 
lem and thereby allow for realistic calculations in man- 
ageable model spaces, we renormalize the interactions 
through a G-matrix procedure for use in the Os-Op-Odls- 
0/lp oscillator basis. Our Hamiltonian is thus given by 
H = t + G((D), where ui is the G-matrix starting energy. 
We use a simple procedure described in Ref. .2JJ] to al- 
leviate the starting-energy dependence of the G-matrix 
in orbitals below the Fermi surface. We also modify the 
Hamiltonian by adding to it the center-of-mass Hamilto- 
nian times a Lagrange multiplier /3 c .m.- Thus, our Hamil- 
tonian becomes H' = H + f3 c . in H CAn ,. We choose /? c . m . 
such that the expectation value of H c , m . — 0.0 MeV. De- 
tails may be found in Ref. |27| . 

We tested the performance of the above coupled cluster 
approximations in the context of the nuclear many-body 
problem by applying them to two closed-shell nuclei, 4 He 
and 16 0, in the one-particle space of four major oscilla- 
tor shells. Shell model diagonalization provided an exact 
answer for a given Hamiltonian in the He case. Compar- 
ing the exact ground- and excited-state energies resulting 
from the diagonalization of the Hamiltonian in the small 
model space with the coupled cluster energies obtained 
in the same model space, we can assess the usefulness 
of various coupled cluster approximations in calculations 
for atomic nuclei. In particular, we can learn about the 
possible role of triply excited clusters in an accurate de- 
scription of ground and excited states without confusing 
the inaccuracies resulting from the inadequate treatment 
of the many-body problem by a given coupled cluster ap- 
proximation with other sources of error. 

We report our results for the ground-state energy of 
4 He in Table [I] We used two types of reference determi- 
nants |$): one constructed from the lowest-energy oscil- 
lator states and the Hartree-Fock determinant. Through- 
out the table we see that the results obtained in the os- 
cillator basis are lower in energy when compared to those 
obtained in the Hartree-Fock basis. The two best meth- 
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ods in the oscillator basis are CR-CCSD(T),a/Ao = 1 
and CR-CCSD(T),c/A = 1. They yield results within 
40 keV and 300 keV of the full configuration interac- 
tion (CI) diagonalization problem, respectively. The CR- 
CCSD(T),a/A = 1 approach applied to the oscillator 
basis overshoots the exact result, which is a consequence 
of using the standard MBPT denominators in the def- 
inition of T 3 , Eq. ©. The CR-CCSD(T),a approach, 
in which the triples correction So is renormalized via the 
presence of the Ao denominator in Eq. (J2J), is more sta- 
ble in this regard, providing the upper bound to the en- 
ergy, although the 600 keV error obtained with the CR- 
CCSD(T),a method is not as impressive as the 40 keV 
error obtained with CR-CCSD(T),a/Ao = 1. In general, 
the CR-CCSD(T) results are considerably more accurate 
than the results of the CCSD calculations, in which T3 is 
ignored, although the CCSD approach describes the bulk 
of the correlation effects, reducing the large 16.273 MeV 
error obtained by calculating with the oscillator 

reference |$) to 1.5 MeV. The effectiveness of the CCSD 
approach can also be illustrated by comparing the CCSD 
energy with the results of truncated shell-model calcula- 
tions (CISD), in which the Hamiltonian is diagonalized in 
the same space of singly and doubly excited determinants 
as used in the CCSD caculations. The costs of the CISD 
and CCSD calculations are almost identical (both are 
v? n\ procedures), and yet the error in the CISD energy 
obtained in the oscillator basis is twice as large as the er- 
ror obtained with CCSD. The noniterative triples correc- 
tions defining the CR-CCSD(T) approaches reduce these 
errors to as little as 40 keV, which is a lot better than 
the 1.3 MeV error in the CISDT calculations, where the 
Hamiltonian is diagonalized in the much larger space of 
all singly, doubly, and triply excited determinants. This 
demonstrates the advantages of coupled cluster methods 
over the diagonalization techniques. Similar observations 
apply to the Hartree-Fock basis, although the coupled 
cluster results obtained with this basis are not as good 
as those obtained with the oscillator basis. For example, 
the best result in the Hartree-Fock basis, obtained with 
CR-CCSD(T),c/A = 1, is 700 keV above the exact re- 
sult. This suggests that we may be better off by using 
the oscillator basis in coupled cluster calculations. On 
the other hand, the CR-CCSD(T) results obtained in the 
Hartree-Fock basis are not unreasonable, allowing us to 
contemplate the use of the Hartree-Fock basis in coupled 
cluster calculations for open-shell nuclei. (This would 
parallel the Hartree-Fock-based coupled cluster calcula- 
tions for open-shell electronic states in chemistry.) 

We used the EOMCCSD method and its CR-CCSD(T) 
extension to compute excited states. To our knowledge, 
this is the first time that nuclear excited states are com- 
puted using coupled cluster methods. The results for 
4 He are given in Table [HI The low-lying J = 1 state 
most likely results from the center-of-mass contamination 
which we have removed only from the ground state. The 



TABLE I: The ground-state energies of 4 He calculated using 
the oscillator (Osc) and Hartree-Fock (HF) basis states. Units 
are MeV. The reference energies ($|//'|<1>) are -7.211 (Osc) 
and -10.520 (HF) MeV. 



Method 




Osc 


HF 


CCSD 




-21.978 


-21.385 


CR-CCSD(T),a 




-22.841 


-22.395 


CR-CCSD(T),a/A = 


= 1 


-23.524 


-22.711 


CR-CCSD(T),b 




-22.396 


-22.179 


CR-CCSD(T),b/A = 


= 1 


-22.730 


-22.428 


CR-CCSD(T),c 




-22.630 


-22.450 


CR-CCSD(T),c/A = 


= 1 


-23.149 


-22.783 


CISD 




-20.175 


-20.801 


CISDT 




-22.235 




Exact 




-23.484 


-23.484 



TABLE II: The excitation energies of 4 He calculated using 
the oscillator basis states (in MeV). 



State EOMCCSD CR-CCSD(T) a CISD Exact 



J= 


1 


11.791 


12.044 


17.515 


11.465 


J= 





21.203 


21.489 


24.969 


21.569 


J= 


2 


22.435 


22.650 


24.966 


22.697 



a The difference of the CR-CCSD(T),c energy of excited 
state and the CR-CCSD(T),b energy of the ground state. 

J = and J = 2 states calculated using EOMCCSD and 
CR-CCSD(T) are in excellent agreement with the exact 
results. For these two states, the EOMCCSD approach 
provides the relatively small, 0.3-0.4 MeV, errors, which 
are further reduced by the CR-CCSD(T) triples correc- 
tions to < 0.1 MeV. Based on the experience with the 
equation of motion coupled cluster methods in chemistry 
|1 II [3 13 EJ S3, the very good performance of the 
EOMCCSD approach for the lowest-energy excited states 
of 4 He can be understood if we realize that these states 
are dominated by single-particle excitations. Again, a 
comparison of the EOMCCSD and CISD results shows 
that coupled cluster theory offers much higher accuracies 
compared to truncated diagonalization of similar numer- 
ical effort. According to the experiment, the lowest ly- 
ing + state in 4 He is a resonance at an excitation of 
20.21 MeV and a width of 0.5 MeV, while the first J = 2 
state lies at 21.84 MeV and has a width of 2 MeV. We 
have not identified the parity of our calculated states, 
but it seems to us that we will be able to model the 
excitations in 4 He and other nuclei using coupled clus- 
ter theory. Our results in Table [H] are indicative of the 
accuracies we may expect from such calculations. 

We also applied the CCSD and CR-CCSD(T) methods 
to 16 O. Table IIIII shows the total ground-state energy 
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TABLE III: The ground-state energy of 16 calculated using 
various coupled cluster methods and oscillator basis states. 
The reference energy {$\H'\$) is -109.452 MeV. 



Method 




Energy 


CCSD 




-139.310 


CR-CCSD(T),a 




-139.465 


CR-CCSD(T),a/A = 


1 


-139.621 


CR-CCSD(T),b 




-139.375 


CR-CCSD(T),b/A = 


1 


-139.440 


CR-CCSD(T),c 




-139.391 


CR-CCSD(T),c/A = 


1 


-139.467 


CISD 




-131.887 


CISDT 




-135.489 


CISDTQ 




-138.387 



values obtained with the CCSD and CR-CCSD(T) ap- 
proaches. As in the 4 He case, coupled cluster methods re- 
cover the bulk of the correlation effects, producing the re- 
sults of the CISDTQ, or better, quality. CISDTQ stands 
for the expensive shell-model diagonalization in a huge 
space spanned by the reference and all singly, doubly, 
triply, and quadruply excited determinants (the most ex- 
pensive steps of CISDTQ scale as n*n^). To understand 
this result, we note that the CCSD T\ and amplitudes 
are similar in order of magnitude. (For an oscillator ba- 
sis, both T\ and T2 contribute to the first-order MBPT 
wave function.) Thus, the T1T2 disconnected triples are 
large, much larger than the T3 connected triples, and the 
difference between the CISDT (CI singles, doubles, and 
triples) and CISD energies is mostly due to T\T-2. The 
small T3 effects, as estimated by CR-CCSD(T), are con- 
sistent with the CI diagonalization calculations. If the T3 
corrections were large, we would observe a significant low- 
ering of the CCSD energy, far below the CISDTQ result. 
The CISDTQ diagonalization is not size-extensive, while 
the CCSD and CR-CCSD(T)/A = 1 methods maintain 
this property. Moreover, the CCSD and CR-CCSD(T) 
methods bring the nonncgligiblc highcr-than-quadruplc 
excitations, such as X1T2, TiXf, and T%, which are not 
present in CISDTQ. It is, therefore, quite likely that the 
CR-CCSD(T) results are very close to the results of the 
exact diagonalization, which cannot be performed. 

In summary, we used the quantum chemistry inspired 
coupled cluster approximations to calculate the ground 
and excited states of the 4 He and 16 nuclei. By com- 
paring coupled cluster results with the exact results ob- 
tained by diagonalizing the Hamiltonian in the same 
model space, we demonstrated that relatively inexpen- 
sive coupled cluster approximations recover the bulk of 
the nuclcon correlation effects in ground- and excited- 
state nuclei. These results are a strong motivation to 
further develop coupled cluster methods for the nuclear 



many-body problem, so that accurate ab initio calcula- 
tions for small- and medium-size nuclei become as routine 
as the molecular electronic structure calculations. 
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